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Abstract 

A  numerical  simulation  tool  for  calculating  the  planar  and  mono-block  layer  built  (MOLB)  type  solid  oxide  fuel  cells  (SOFC)  is  described. 
The  tool  combines  the  commercial  computational  fluid  dynamics  simulation  code  with  an  electrochemical  calculation  subroutine.  Its  function  is  to 
simulate  the  heat  and  mass  transfer  and  to  predict  the  temperature  distribution  and  mass  fraction  of  gaseous  species  in  the  SOFC  system.  The  three- 
dimensional  geometry  model  of  SOFC  was  designed  to  simulate  a  co-flow  case  and  counter-flow  case.  The  finite  volume  method  was  employed 
to  calculate  the  conservation  equations  of  mass,  momentum  and  energy.  Moreover,  the  influences  of  working  conditions  on  the  performances  of 
planar  and  MOLB-type  SOFCs  were  also  discussed  and  compared,  such  as  the  delivery  rate  of  gas  and  the  components  of  fuel  gas.  Simulation 
results  show  that  the  MOLB-type  SOFC  has  higher  fuel  utilization  than  the  planar  SOFC.  For  the  co-flow  case,  average  temperatures  of  PEN 
(positive  electrode-electrolyte-negative  electrode)  in  both  types  of  SOFCs  rise  with  the  increase  in  delivery  rate  and  mass  fraction  of  hydrogen.  In 
particular,  the  temperature  of  planar  SOFC  is  more  sensitive  to  the  working  conditions.  In  order  to  decrease  the  average  temperatures  in  SOFC,  it 
is  effective  to  increase  the  delivery  rate  of  air. 

©  2007  Elsevier  B.V.  All  rights  reserved. 

Keywords:  Planar  SOFC;  MOLB-type  SOFC;  Thermo-electrochemical  model;  Temperature  distribution;  Gas  mass  fraction  distribution 


1.  Introduction 

Solid  oxide  fuel  cell  has  been  looked  as  a  promising  alterna¬ 
tive  energy  source  for  residential  and  distributed  power  plants 
because  of  its  higher  energy  conversion  efficiency  and  power 
density,  low  environmental  hazards  and  potentially  low  produc¬ 
tion  cost.  Thus,  SOFC  is  expected  to  realize  commercialization 
within  a  few  years  [1-3].  However,  the  further  development  of 
the  SOFC  has  the  challenges  related  to  maximize  the  power 
density  and  minimize  the  non-uniform  temperature  distribution, 
which  contributed  to  the  thermal  stress  in  the  SOFC  compo¬ 
nents  [4,5],  so  the  heat  and  mass  transfer  in  the  SOFC  need  to 
be  researched. 
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The  heat  and  mass  transfer  in  the  SOFC  are  influenced  not 
only  by  geometrical  designs  of  SOFC  but  also  by  the  work¬ 
ing  conditions.  In  the  past,  some  researchers  have  performed 
on  preparing  the  electrode  and  electrolyte  materials  [6-9]  and 
designing  geometry  of  SOFC  [10-13].  Thereinto,  the  planar 
SOFC  has  received  much  more  attention  than  the  tubular  SOFC 
because  of  its  easier  fabrication  and  higher  power  density.  In 
recent  years,  another  type  SOFC,  which  was  MOLB-type  SOFC, 
has  also  been  researched.  Such  design  has  been  in  favor  of 
enhancing  its  mechanical  strength  and  the  corrugate-shaped 
PEN  provides  the  film  with  the  combined  function  of  fuel  and  air¬ 
flow  paths,  so  the  cell  stacks  are  more  compact  and  the  laborious 
work  of  channel  machining  is  less. 

Besides  the  geometrical  designs,  the  working  conditions, 
such  as  delivery  rates  and  hydrogen  mass  fraction  of  fuel 
and  air  to  the  cell  system  also  influence  on  the  performances 
of  SOFC  in  a  complex  way.  Thus,  in  order  to  efficiently 
develop  SOFC  stacks,  it  is  convenient  and  effective  to  have 
the  capability  to  experiment  numerically  with  the  effects  of 
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geometric  designs  and  operation  conditions  on  SOFC  perfor¬ 
mances.  In  the  past,  modeling  of  the  SOFC  during  steady 
operation  has  been  constructed  to  calculate  temperature  distribu¬ 
tion  [14-20].  Investigations  of  planar  and  MOLB-type  SOFCs 
operation  and  performance  have  predicted  cell  temperature  dis¬ 
tributions  for  various  flow  patterns  [21-24].  However,  only  a 
little  work  has  been  performed  on  the  influences  of  the  oper¬ 
ating  conditions  on  the  performances  of  planar  SOFC  [25,26], 
or  even  saying  nothing  of  the  influences  on  the  MOLB-type 
SOFC. 

The  objective  of  present  work  is  to  compare  the  heat  and  mass 
transfer  in  the  conventional  planar  with  that  in  the  MOLB-type 
SOFC.  A  CFD  model  tool  is  demonstrated  to  predict  temperature 
distributions  and  gas  species  distributions  of  two-design  SOFCs 
system.  The  results  simulated  in  this  paper  can  not  only  guide  the 
designer  in  understanding  how  geometrical  design  affects  the 
thermodynamics  performances  in  the  SOFC,  but  also  provide 
a  more  reasonable  basis  for  geometrical  design  of  the  SOFC 
stacks. 

2.  Mathematics  model 

2.7.  Model  geometry 

Fig.  1(a)  schematically  shows  a  typical  planar  SOFC.  PEN 
plate  is  sandwiched  between  the  inter-connectors,  which  are 
machined  with  passages  for  channeling  the  fuel  and  air  to  the 


fuel  cell.  Fig.  1(b)  demonstrates  the  MOLB-type  SOFC  and  the 
PEN  is  modeled  into  a  corrugated  shape  [21].  The  fuel  and  air 
flow  in  the  trapezoidal  channels  are  separated  by  the  corrugated 
PEN.  For  sake  of  simplicity  in  the  calculation,  one  repeating 
cell  unit  was  analyzed  in  this  simulation  (as  shown  in  Fig.  1).  In 
this  model,  for  the  two-design  SOFCs,  the  thicknesses  of  anode, 
cathode,  electrolyte  and  inter-connector  were  0.5, 0.25, 0.05  and 
1 .0  mm,  respectively. 

2.2.  Thermo -fluid  model 

The  ANSYS-CFX  code  was  selected  to  solve  the  thermo¬ 
fluid  model.  In  the  simulation,  the  conservation  equations  of 
species,  mass,  momentum  and  energy  were  solved  using  the 
finite  volume  method. 

In  general,  gas  species  transfer  mainly  by  convection  in  the 
flow  channels  and  diffusion  in  the  porous  electrodes.  The  species 
conservation  equation: 

V(pCkU)  =  V(D*feffVC*)  +  Ik,  k  =  H2,  02,  H20  (1) 

where  4  is  the  rate  of  production  or  consumption  of  species  k, 
and  given  by  [27] : 

4  =  ±g;  (2) 

where  Sk  is  the  chemical  computation  coefficients,  i  is  the  local 
current  density,  and  F  is  the  Faraday  constant. 


Fig.  1.  Illustrations  of  the  one  cell-stack  and  the  single  unit  model  for  planar  SOFC  (a)  and  MOLB-type  SOFC  (b). 
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The  diffusion  coefficient  of  gas  is  obtained  by  [28]: 

/  T  V'5 

£>h2  =  0.753  x  lOe  -4  x  (  —  j  , 

Dq2  =0.181  x  10e-4x  (3) 

where  Dh2  and  Do2  are  the  diffusivity  of  H2  and  O2,  respec¬ 
tively,  unit  is  m2  s_1. 

Although  the  mass  is  added  on  the  anode  and  removed  on  the 
cathode,  but  the  decrement  of  mass  on  the  cathode  is  equal  to 
increment  of  mass  on  the  anode,  so  the  mass  in  the  whole  SOFC 
system  is  conservational.  The  mass  conservation  equation: 


X(spU)  =  0  (4) 

Both  the  air  and  fuel  flows  were  considered  as  ideal  gas  mixtures 
with  the  density  given  by 


Table  1 


Coefficients  of  the  specific  heats  of  gas 


Gas 

a 

b 

c 

Hydrogen 

25.8911 

-0.8373 

2.0138 

Oxide 

29.0856 

12.9874 

-3.8644 

Water  gas 

30.3794 

9.6212 

1.1848 

Heat  transfer  between  the  fluid  and  solid  materials  was  lim¬ 
ited  to  conduction  and  convection,  and  radiation  was  neglected 
in  this  calculation  because  it  is  very  small  relative  to  the  other 
kinds  of  heat  transfer  mechanisms.  So  the  energy  conservation 
equation  is: 

V(t/(p£f  +  p ))  +  V(r  U)  +  V(keffVT)  +  SE  =  0  (10) 

where  p  is  the  density  of  the  gas  fluid,  Ef  is  the  total  fluid  energy, 
t  is  the  stress  tensor,  and  keff  is  the  effective  thermal  conduc¬ 
tivities  of  porous  electrodes,  which  are  calculated  as  follows 
[29]: 


where  mg  is  the  mass  fraction  of  species  k  with  molecular  weight 
Mg  and  R  is  the  gas  constant. 

The  momentum  conservation  equation: 


ps 


3  u  dv  3  w 

u - b  v - b  w  — 

dx  3  y  3  z 


=  -sVP  +  £/Xe  ff 


/  3  u2  dv 2  3  w2\ 

\  dx2  3y2  3z2  ) 


+  Sm 


(6) 


where  s  is  the  porosity  of  electrode,  P  is  the  pressure,  /xeff  is  the 
effective  viscosity  of  the  mixture  gas  and  is  given  by  [27] : 


/4eff  =  Y 
k 


XkP'k 

'f2Xk<t>ki 

k 


(7) 


.  [1  +Wk/l^j)l/\MjIMk)l/4]2 

<Pki  = - - - rjTj -  (») 

[8(1  +  Mk/Mj)] 1/2 

In  Eq.  (6),  Sm  is  momentum  source  and  Sm  =  0  in  the  flow  chan¬ 
nels.  However,  in  the  porous  electrodes,  Darcy  law  with  constant 
porosity  and  permeability  is  applied  to  model  the  momentum 
source  as  follows  [25]: 


Su  =  -^s2U  (9) 

K 

where  Xg  is  the  molar  percent  of  the  species  k ,  /xy,  and  pig  are 
kinematical  viscosities  of  species  j  and  k ,  respectively. 


fceff  =  ski  +  (1  -  e)ks  (11) 

where  kf  and  ks  are  thermal  conductivities  of  fluid  and  solid, 
respectively.  Se  is  energy  source  expressed  by  Eq.  (12)  and 
mainly  consists  of  reaction  and  Ohmic  heats  [27] . 

i2  i  ( TAS  \ 

Se  — - b  -  (  +  rj  act  )  (12) 

°eff  8  \  2F  J 

where  i  is  the  local  current  density,  creff  is  the  effective  electri¬ 
cal  conductivity,  8  is  the  anode  thickness,  AS  is  the  change  of 
entropy  in  the  reaction  and  ?7act  is  the  activation  potential. 

Temperature-dependent  specific  heat  of  gas  is  used  in  the 
simulation,  which  are  as  follows  [27]: 

Cp  =a  +  b  X  1037’  +  c  x  1067’2  (13) 

where  a ,  b ,  and  c  are  relevant  coefficients,  as  listed  in  Table  1. 
Solid  material  properties  used  in  this  simulation  are  listed  in 
Table  2. 

2.3.  Electrochemical  model 

2.3.1.  Assumptions  and  reactions 

The  oxidant  reduction  reaction  occurring  at  the  cathode  is 
expressed  as  follows: 

0.5  02  +  2e“  -*  O2-  (14) 


Table  2 

Properties  of  the  solid  materials 


Cell  component 

Density  (kg  m  3) 

Effect  thermal  conductivity  (W  (m  K)  1 ) 

Specific  heat  (J  (kg  K)-1) 

Porosity  (%) 

Permeability  coefficient  (m2) 

Interconnect 

7700 

13 

0.8 

- 

- 

Anode 

6200 

6.23 

0.65 

35 

1.0E-12 

Cathode 

6000 

9.6 

0.9 

35 

1.0E-12 

Electrolyte 

5560 

2.7 

0.3 

- 

- 
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The  oxygen  ions  transfer  through  the  electrolyte  and  then  into 
the  active  reaction  areas  of  anode.  The  electrochemical  reaction 
of  fuel  at  the  anode  is 

H2  +  O2”  H20  +  2e~  (15) 

So  the  overall  reaction  is 

H2  +  \02  ->  H20  (16) 


2.3.2.  Dynamics  of  electrochemical  reactions 

According  to  the  Faraday  law,  the  reaction  rates  depend  on 
the  current  density  i  [  15]: 


d  /  d02 

i  =  2F—  =4  F — - 
d  t  d  t 


(17) 


where  df/dt  and  d02/df  are  the  molar  consumption  rates  of  fuel 
and  oxygen  at  the  anode  and  the  cathode,  respectively. 

During  the  process  of  energy  transforming,  when  the  charge 
transfer  reaction  at  the  electrolyte-electrode  interface  is  too 
slow  to  provide  ions  at  the  rate  required  by  the  demand  of  cur¬ 
rent,  the  activation  polarization  occurs  and  is  defined  by  the 
Butler- Volmer  equation  [30] : 


i  =  i0 


IF 

exp(l  -  A)  —  /? act 


(18) 


Theoretical  voltage  of  cell  is  obtained  by  the  Nemst  equation: 


RT 

E0  H - In 

U  IF 


(  Ph2o  \ 
V Pu2  x  Pol) 


RT 


= - In  K  +  RT  x  In 

IF 


PHl  X  Pol \ 

Pu2o  J 


(21) 


where  K  is  the  equilibrium  constant  of  reaction  under  the  stan¬ 
dard  atmosphere,  Eq  is  the  standard  voltage  of  the  cell,  /?h2o, 
po2 ,  and  pn2  are  the  partial  pressures  of  water  gas,  oxygen  and 
hydrogen,  respectively. 


3.  Numerical  implementation 

In  the  calculations,  the  modeling  tool  couples  an  electro¬ 
chemical  calculation  method  with  a  commercial  computational 
fluid  dynamics  (CFD)  simulation  code.  The  finite  volume 
Navier-Stokes  and  transport  equations  are  solved  to  obtain  the 
gas  species  mass  fraction  and  temperatures  at  each  position  in 
the  cell.  The  information  is  passed  to  the  electrochemical  model 
(subroutine).  Then  the  local  current  density  is  calculated  and 
applied  to  obtain  the  hydrogen  reaction  rate,  heat  source  and 
species  sources.  Gas  species  mass  fraction  and  temperature  dis¬ 
tributions  are  then  calculated  for  the  next  iteration,  and  so  on, 
until  convergence  of  solution  is  achieved. 


Eq.  (18)  is  simplified  and  described  by  the  empirical  formula 
and  Tafel  equation  in  anode  and  cathode,  respectively  [29,30]: 


4.  Simulation  results  and  discussion 


7? act,  a  — 


RT 

2^0,a 


(19) 


/  RT  \  (  RT  \ 

•**u=-{w)'r‘i0-c+{w)'"i  m 

where  f  is  the  transmission  coefficient  and  ft  =  0.5  in  this  simula¬ 
tion,  Osict, a  and  ?7act,c  are  the  activation  potentials  at  the  anode  and 
the  cathode,  respectively.  z'o,a  and  z'o,c  are  the  exchange  current 
densities  at  the  anode  and  the  cathode,  respectively. 


The  cell  voltage  imposed  on  both  geometries  is  0.3  V.  In  the 
porous  electrodes,  typical  values  are  employed  for  porosity  of 
both  electrodes  £  =  0.5,  tortuosity  of  anode  is  4.5,  tortuosity  of 
cathode  is  3,  and  permeability  is  1.7  x  10_10m2.  In  this  study, 
exchange  current  densities  at  the  anode  and  the  cathode  are 
5300  A  m-2  and  2300  A  m-2,  respectively.  And  other  cell  oper¬ 
ating  conditions,  including  the  delivery  rate  of  fuel  and  hydrogen 
mass  fraction  in  the  fuel  inlet  are  listed  in  Table  3. 

Fig.  2  compares  the  PEN  temperature  distributions  between 
planar  SOFC  and  the  MOLB-type  SOFC  for  co-flow  case  (no.  1 ). 


Table  3 

The  cell  operating  conditions  and  parameters  used  for  simulation 


Sample  number  case 

Fuel 

Air 

Flow  pattern 

Delivery  rate 

Inlet  temperature  (K) 

Hydrogen  mass 

Delivery  rate 

Inlet  temperature  (K) 

(vl)  (ms-1) 

fraction  (%) 

(v2)  (ms-1) 

1 

0.5 

973 

0.8 

3 

873 

Counter-flow 

2 

0.5 

Co-flow 

3 

1.0 

973 

0.8 

3 

873 

Co-flow 

1.5 

0.8 

4 

0.5 

973 

0.9 

1 

3 

873 

Co-flow 

3 

5 

0.5 

973 

0.8 

2 

1 

873 

Co-flow 
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Temperature  [K] 
1.053e+003 

1 . 038e+003 

1 . 023e4-003 

1.00764-003 

9. 923e*t-002 

9. 771e4-002 

9. 6 19e4002 

9.467e4-002 

9. 315e4-002 

9. 163e4002 


(b) 


Temperature  [K] 
-  1  063C4003 


-  1  054C4003 

-  1.045e+003 

-  1.03764003 
1.028e4003 

L  I.019e4003 
L  1.010e4003 


1 . 00 1e4003 
9  925e4002 
9  837e4002 


Fig.  2.  PEN  temperature  distributions  for  planar  SOFC  (a)  and  MOLB-type  SOFC  (b)  in  the  co-flow  case  (no.  1). 


Although  the  structure  of  SOFC  between  these  two  geometries  is 
different,  the  general  trends  in  the  temperature  distributions  are 
similar  in  essential.  It  is  briefly  described  as  follows.  The  aver¬ 
age  PEN  temperature  in  planar  SOFC  is  979  °C  with  maximum 
and  minimum  temperatures  of  1053  °C  and  916  °C,  respectively. 
Moreover,  the  PEN  average  temperature  is  1014  °C  with  max¬ 
imum  and  minimum  temperatures  of  1063  °C  and  987  °C  in 
MOLB-type  SOFC.  Thereby,  the  average  temperature  is  higher 
and  the  temperature  difference  (AT)  is  lower  in  MOLB-type 
SOLC.  This  is  because  the  electrochemical  reaction  active  area 
of  MOLB-type  SOLC  is  larger  than  that  of  planar  SOLC,  so  more 
gas  is  consumed  in  the  active  areas,  as  a  result,  more  reaction 
heat  is  accumulated  and  the  average  temperature  is  higher. 

Lig.  3  compares  the  PEN  hydrogen  mass  fraction  distribu¬ 
tions  between  planar  SOLC  and  the  MOLB-type  SOLC  for 
co-flow  case  (no.  1).  The  hydrogen  mass  fraction  for  both  types 
of  SOLCs  decrease  along  the  fuel  flow  direction,  but  the  dif¬ 
ference  of  hydrogen  mass  fraction  in  the  fuel  inlet  and  outlet 
is  lower  in  MOLB-type  SOLC  (Lig.  3(b))  than  that  in  planar 
SOLC  (Lig.  3(a)),  So  more  hydrogen  is  consumed  because  the 
hydrogen  mass  fraction  in  the  fuel  inlet  is  same,  which  means 


that  the  fuel  utilization  is  higher  in  the  MOLB-type  SOLC.  In 
particular,  due  to  the  corner  effect  formed  by  the  inclined  plane 
and  the  upper  plane,  the  temperature  and  hydrogen  mass  fraction 
distributions  in  MOLB-type  SOLC  are  less  uniform  than  those 
in  planar  SOLC. 

Lig.  4  compares  the  temperature  distributions  of  the  two  types 
of  SOLCs  in  the  counter-flow  case  (no. 2),  the  average  temper¬ 
ature  is  996  °C  with  maximum  and  minimum  temperatures  of 
1088  °C  and  923  °C  in  planar  SOLC,  and  the  PEN  average  tem¬ 
perature  is  1014  °C  with  maximum  and  minimum  temperatures 
of  1063  °C  and  987  °C  in  MOLB-type  SOLC.  In  addition,  it 
should  be  noted  that  the  temperatures  of  PEN  rise  rapidly,  reach¬ 
ing  a  maximum  near  the  fuel  inlet,  and  then  gradually  drop  in 
those  two  types  of  SOLCs.  This  is  due  to  offsetting  effects  of  air 
near  the  inlet,  at  its  coolest,  being  aligned  with  the  fuel  inlet. 
Of  the  two  types  of  SOLCs,  the  MOLB-type  SOLC  has  the 
higher  average  temperature  and  lower  temperature  difference 
( AT)  than  that  of  planar  SOLC,  which  is  similar  to  the  co-flow 
case.  Lurthermore,  the  hydrogen  mass  fraction  for  both  types 
of  SOLCs  also  decrease  along  the  fuel  flow  direction,  but  the 
difference  of  hydrogen  mass  fraction  in  the  fuel  inlet  and  outlet 


(a)  Hydrogen  Mass  Fraction 
' 6.800e- 001 

6. 156e-001 

-5.5116-001 

-  4.8676-001 
4.222e-001 

•  3. 5786-001 

*  2. 933e-001 
-2.2896-001 

-  1.644e-001 

-  1.0006-001 


(b) 


Hydrogen  Mass  Fraction 
6. 571e-001 

5 . 896e-001 

5.221e-001 

4.5466-001 

3.8716-001 

3. 195e-001 

2.520e-001 
L  1.845e-001 
-  1.1706-001 


-  4.952C-002 


Fig.  3.  Hydrogen  mass  fraction  distributions  on  the  interface  between  anode/electrolyte  for  the  planar  SOFC  (a)  and  MOLB-type  SOFC  (b)  in  the  co-flow  case 
(no.  1). 
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Fig.  4.  PEN  temperature  distributions  for  planar  SOFC  (a)  and  MOLB-type  SOFC  (b)  in  the  counter-flow  case  (no.  2). 


is  lower  in  MOLB-type  SOFC  (Fig.  5(b))  than  that  in  planar 
SOFC  (Fig.  5(a)),  So  more  hydrogen  is  consumed  because  the 
hydrogen  mass  fraction  in  the  fuel  inlet  is  same,  in  other  words, 
the  fuel  utilization  is  higher  in  the  MOLB-type  SOFC,  which  is 
similar  to  that  in  the  co-flow  case.  However,  the  hydrogen  mass 
fraction  in  the  fuel  outlet  is  higher  in  counter-flow  case  for  two 
types  of  SOFCs,  so  the  less  hydrogen  reacted  and  hydrogen  uti¬ 
lizations  are  lower  than  those  in  co-flow  case  for  the  two  types 
of  SOFCs. 

Consequently,  overall  considering  the  temperature  distribu¬ 
tion  and  the  hydrogen  utilization,  it  is  advantage  and  suitable  to 
choose  the  co-flow  case  for  SOFCs  steady  operating. 

For  this  flow  case,  the  voltage  is  assumed  constant  and  it 
is  0.3  V,  several  parameters  influenced  on  the  performances  of 
two  types  of  SOFCs  are  also  compared.  First,  we  focus  on  the 
effects  of  the  fuel  delivery  rate  on  temperature  distributions  in 
those  types  of  SOFCs.  Fig.  6  shows  the  temperature  distribution 
of  the  mid-plane  in  the  X-direction  for  the  co-flow  case  (no.  3). 
With  the  increase  in  the  delivery  rate  of  fuel,  the  average  tem¬ 
peratures  and  temperature  differences  (AT)  of  PEN  rise,  which 
may  cause  larger  thermal  stress.  This  is  because  increment  of 


the  fuel  delivery  rate  may  cause  the  increasing  of  reaction  rate, 
so  more  reaction  heat  was  accumulated  and  the  cell  temperature 
rise. 

Consistent  with  increasing  of  the  hydrogen  proportion  from 
80%  to  90%  and  100%  in  the  fuel  gas,  the  temperature  differ¬ 
ences  and  the  average  temperatures  also  rise  due  to  the  increase 
of  reaction  rate  and  more  reaction  heat  accumulated,  as  shown 
in  Fig.  7.  Fig.  7  shows  the  temperature  distribution  of  the  mid¬ 
plane  in  the  X-direction  for  the  co-flow  case  (no.  4),  although  the 
general  variation  trends  of  temperatures  are  similar,  the  temper¬ 
ature  of  planar  SOFC  is  more  sensitive  to  the  working  conditions 
due  to  the  less  electrochemical  reaction  areas. 

On  the  basis  of  the  analysis  on  the  temperature  distributions 
mentioned  above,  in  order  to  decrease  the  average  temperatures 
and  maximum  temperatures,  it  is  effective  to  increase  delivery 
rate  of  air.  In  the  SOFC  system,  air  not  only  provides  oxygen  ions 
but  also  has  the  cooling  function.  Increasing  the  delivery  rate  of 
air,  the  average  and  maximum  temperatures  in  those  types  of 
SOFCs  drop,  as  shown  in  Fig.  8.  This  is  because  more  reaction 
heat  is  absorbed  and  released  by  the  air  with  the  higher  delivery 
rate,  although  the  air  utilizations  are  dropped. 
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Fig.  5.  Hydrogen  mass  fraction  distributions  on  the  interface  between  anode/electrolyte  for  the  planar  SOFC  (a)  and  MOLB-type  SOFC  (b)  in  the  counter-flow  case 
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Fig.  6.  Comparison  of  PEN  temperature  distributions  between  planar  SOFC  and  MOLB-type  SOFC  in  the  co-flow  case  (no.  3). 
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Fig.  8.  Comparison  of  PEN  temperature  distributions  between  planar  SOFC  and  MOFB-type  SOFC  in  the  co-flow  case  (no.  5). 
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5.  Conclusions 

The  local  current  density  was  used  to  couple  the  thermo-fluid 
model  with  electrochemical  model.  The  temperature  distribu¬ 
tions  of  planar  and  MOLB-type  SOFCs  in  steady  state  were 
compared  and  some  main  conclusions  were  made  as  follows: 

(1)  In  the  co-flow  case,  the  PEN  temperature  distributions  in 
planar  and  MOLB-type  SOFCs  are  more  uniform  than  those 
in  counter-flow  case. 

(2)  Under  the  same  work  conditions,  the  temperature  gradient  is 
lower  in  MOLB-type  SOFC  than  those  in  the  planar  SOFC. 

(3)  For  co-flow  case,  under  the  constant  voltage  condition,  with 
increasing  the  delivery  rate  of  fuel  gas  or  hydrogen  mass 
fraction  in  the  fuel,  temperature  gradients  in  the  two  types 
of  SOFCs  rise.  Moreover,  the  temperature  distributions  in 
the  planar  SOFC  are  more  sensitive  to  the  variation  of  deliv¬ 
ery  rate  of  fuel  and  hydrogen  mass  fraction.  In  an  attempt 
to  achieve  a  more  uniform  temperature  distribution,  it  is 
effective  to  decrease  the  temperature  gradients  of  PEN  by 
increasing  the  delivery  rate  of  air. 
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